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We consider a one-dimensional system consisting of electrons with short-ranged repulsive inter- 
actions and coupled to small-momentum transfer acoustic phonons. The interacting electrons are 
bosonized and described in terms of a Luttinger liquid which allows to calculate exactly the one- 
and two-electron Green function. For non-interacting electrons, the coupling to phonons alone in- 
duces a singularity at the Fermi surface which is analogous to that encountered for electrons with an 
instantaneous attractive interaction. The exponents which determine the presence of singlet /triplet 
superconducting pairing fluctuations, and spin/charge density wave fluctuations are strongly affected 
by the presence of the Wentzel-Bardeen singularity, resulting in the favoring of superconducting fluc- 
tuations. For the Hubbard model the equivalent of a phase diagram is established, as a function of: 
the electron-phonon coupling, the electron fllling factor, and the on site repulsion between electrons. 
The Wentzel-Bardeen singularity can be reached for arbitrary values of the electron-phonon cou- 
pling constant by varying the filling factor. This provides an effective mechanism to push the system 
from the antiferromagnetic into the metallic phase, and finally into the superconducting phase as 
the electron filling factor is increased towards half filling. 
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I. INTRODUCTION 



In the light of recent technological advances in nanos- 
tructure fabrication, there has been a renewed interest 
in the properties of one-dimensional electron systems. A 
two-dimensional electron gas created in a GaAS / AlGaAs 
heterostructure can be confined in its lateral dimension 
by means of metallic gates, yielding to a good approxima- 
tion a one-dimensional wire. Up to recently, the theoret- 
ical study of quasi-one-dimensional mesoscopic systems 
has focused mostly on non-interacting electrons. In 
the present work, we will consider the properties of a sys- 
tem of electrons in one dimension interacting via a short 
ranged repulsive potential, which in turn interacts with 
phonons through a deformation potential coupling. The 
coupling to the phonons gives rise to a retarded attractive 
electron-electron interaction which opens the possibility 
for superconducting fluctuations. On the other hand, 
the instantaneous repulsion between electrons may be so 
strong that only spin density wave (SDW) and charge 
density wave (CDW) fluctuations can survive. In addi- 
tion, the type of fluctuations which dominate will depend 
crucially on the electron filling factor. We will show that 
the interplay between the phonon-mediated attractive 
interaction and the instantaneous repulsive interaction 
can be studied systematically in one dimension, and thus 
obtaining the equivalent of a phase diagram for the fluc- 
tuations. 

Aside from the experimental motivations mentioned 
above, one-dimensional models can provide useful an- 
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swers: they are in general easier to handle than their 
higher dimensional analogues, for which a full treatment 
of the correlation effects is often difficult. While the one- 
dimensionality strictly rules out the possibility of long 
range order, it is still legitimate to ask what type of or- 
dering fluctuations exist in such systems, such as singlet 
and triplet pairing fluctuations, or spin density wave and 
charge density wave fluctuations. This question was ad- 
dressed for a purely instantaneous interaction between 
electrons in the seventies and later on for non- 

interacting electrons coupled to phonons ||^. More re- 
cently, the attention has focused on interacting electrons 
coupled to phonons |^,^J|J^. 

The interest in the latter systems is by no means new: 
static properties of ID interacting electron systems cou- 
pled to low energy acoustic phonons were considered four 
decades ago by Wentzel |0 and Bardeen jl^] as a pos- 
sible candidate for the theory of superconductivity. The 
electrons were described in terms of Tomonaga bosons 
| p2[ , and it was shown that for a critical value of the 
electron-phonon coupling constant, the system becomes 
unstable and aquires a negative compressibility. The 
thermodynamic quantities for the electron-phonon sys- 
tem (within the Einstein model) were studied near this 
singular point |22j, indicating the presence of a phase 
transition: the divergence of the specific heat is accom- 
panied by a collapse of the system induced by the strong 
electron-phonon interaction. We shall refer to this point 
as the Wentzel-Bardeen singularity in the course of this 
work. Nevertheless, the Wentzel-Bardeen singularity was 
not taken too seriously in the context of superconductiv- 
ity, as the (strong) coupling needed to trigger this dra- 
matic behavior was judged to have unphysical values. 

More recently, a system similar to that of ref. |2|] was 
examined, with the conclusion that the coupling to low 
energy acoustic phonons can be neglected in most sys- 
tems because the effects associated with the phonons con- 
tain the small parameter c/vp, where c is the speed of 
sound and vp is the Fermi velocity. As a result, this 
type of coupling ceased to receive much attention, and 
the interest shifted towards higher momentum processes 
where phonons transfer (or backscatter) an electron from 
one side of the Fermi surface to the other side 0,^,^: 2kp 
phonon processes. These processes are understood to be 
at the origin of the Peierls instability. Nevertheless, in the 
present work, we shall come back to the case of low energy 
acoustic phonons, and will point out that for strongly 
correlated electron systems, this type of coupling can 
have spectacular effects as one approaches the Wenztel- 
Bardeen singularity. For strongly correlated systems, the 
effects associated with the phonons do not depend on the 
small parameter c/vp: rather, the Fermi velocity vp is 
replaced by the charge velocity Up associated with low en- 
ergy particle-hole excitations close to the Fermi surface. 
In contrast to vp, Up can depend strongly on the electron 
filling factor: for the Hubbard model, Up vanishes at half 
filling whereas vp does not. Near half filling, the decay 
of Up is most dramatic for low values of the on site repul- 
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sion parameter U, i.e. the system is extremely sensitive 
to the fiUing factor. As a consequence, the effects associ- 
ated with these low energy phonons can no longer be dis- 
missed as a small correction. Below, we will argue that: 
1) the Wentzel-Bardeen singularity can be reached for 
arbitrary non-zero value of the electron-phonon coupling 
constant as one approaches half filling in the Hubbard 
model. 2) Near the Wentzel-Bardeen singularity, CDW 
and SDW fluctuations are strongly suppressed, and the 
system is first pushed towards an intermediate (metallic) 
phase, and eventually into the superconducting phase. 

The low energy properties of a one-dimensional elec- 
tron system can be treated rigorously using the so-called 
bosonization technique. This approach has been applied 
recently to a variety of problems in mesoscopic physics: 
the transport properties of quasi-one- dimensional inter- 
acting electrons systems ||l^ , [l7|] , the persistent current 
of a one-dimensional ring of interacting fermions [ p^Jl9| 
among others. In the Luttingcr liquid picture, the field 
operator describing the electrons are characterized by a 
charge field and a spin field, which refer to the two differ- 
ent types of elementary excitations which are present in 
the system at low temperatures. In the Luttinger liquid 
approximation, these excitations are totally decoupled, 
and each is characterized by a separate velocity. Because 
of this, charge and spin excitations can be localized at 
different points in space. For the particular type of in- 
teraction considered here, the coupling to the phonons 
occurs only with the charge degrees of freedom. The 
collective excitations describing the charge-phonon dy- 
namics are then characterized by two distinct velocities, 
which depend on the charge velocity, the sound veloc- 
ity, and the electron-phonon coupling. This property is 
reminiscent of the phenomenon of hybridization. As a 
consequence, the electronic system coupled to phonons is 
characterized by three distinct velocities. 

The first part of this paper will be concerned with the 
calculation of the single-electron Green function where 
the phonon part is integrated out. The correspond- 
ing Green function for interacting spinless fermions cou- 
pled to phonons has been calculated in the early eight- 
ies within perturbation theory, assuming a linear spec- 
trum near the Fermi surface |23|| , and the results were ex- 
panded in the ratio c/vp- Here, we shall use a functional 
integral approach, where the phonons can be explicitly 
traced out of the problem. This is the equivalent of sum- 
ming up all phonon diagrams of the perturbation series 
p3| . The results for the Green function can in turn be 
used to determine the the momentum distribution func- 
tion of the electrons, which is of the Fermi-Dirac type 
for a free electron gas. However, both the instantaneous 
interaction and the phonon-mediated retarded interac- 
tion contribute to drive the system away from the free 
electron behavior. This is in analogy with the early work 
on bosonization which showed that for one-dimensional 
electron systems, the momentum distribution function 
departs from the Fermi-Dirac distribution for arbitrarily 
small coupling p4| . 
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In a second part, we will focus on the many body prop- 
erties of the system. Two-particle correlation functions 
in Matsubara representation can be computed with the 
same method, yielding the long range and long (imag- 
inary) time behavior of the singlet and triplet pairing 
correlation function, as well as of the the SDW and 
CDW correlation functions. For an uncorrelated fermion 
system in the limit of vanishing temperature, arbitrar- 
ily small phonon coupling leads to a divergence of the 
Fourier transform of the pairing correlation function at 
low momentum and low frequency, which is precisely the 
signature for superconducting fluctuations This pic- 
ture will be modified when we include repulsive interac- 
tions. For an interacting electron system with a short 
range repulsive potential between electrons and with a 
filling factor between zero and half filling, but not close 
to these values, the onset of superconducting fluctuations 
correspond in general to a finite value of the electron- 
phonon coupling constant. This threshold value of the 
phonon coupling is directly related to the parameters of 
the electronic model which is considered. As an example, 
we study the Hubbard model, where the Luttinger liquid 
parameters can be computed exactly. This will allow us 
to study the instabilities of the system for filling factors 
close to half filling, and therefore to study the sensibility 
to the Wentzel-Bardeen singularity. 

Finally, we can apply our results to study the behavior 
of two coupled electron system or alternatively a electron 
system coupled to a hole system, with a local interaction. 
The density-density coupling has then precisely the same 
form as for electrons coupled to low energy phonons, and 
consequently there is the possibility that one system may 
induce superconducting fluctuations in the other. While 
it is possible that the inclusion of 2kF phonons could 
modify the behavior of the elect ron-phonon system, for 
two coupled electron systems, there is no such possibil- 
ity: the Fermi velocities of each mode arc fixed and dis- 
tinct from each other, and there is no cancellation of the 
fast varying components of the density operator. Only 
the slow varying component of these operators has to 
be taken into account, and our model is "exact" in this 
sense. 

The paper is organized as follows: in Sec. II, the 
Hamiltonian of the coupled electron-phonon system is 
described, and a brief review of the bosonization tech- 
nique is presented. In Sec. Ill, we introduce the func- 
tional integral formalism to calculate the thermodynamic 
quantities. Sec. IV is devoted to the calculation the 
single particle Green function. The ordering correlation 
functions are considered in Sec. V. We illustrate our re- 
sults with the Hubbard model in Sec. VI. The relevance 
of our results for coupled electron chains are discussed in 
Sec. VII. We conclude in Sec. VIII. 
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II. DESCRIPTION OF THE MODEL 



The starting point of our approach is a system of inter- 
acting electrons on a lattice coupled to acoustic phonons 
in one dimension, with periodic boundary conditions. 
The Hamiltonian describing the acoustic phonons is given 
by: 

Hp^^Jdx [C'Ul + Cc'id.d)'] , (2.1) 

where d is the phonon field operator, is its canonical 
conjugate, C is the mass density of the lattice and c is the 
speed of sound. 

Let "08 be the field operator for fermions with spin s. In 
the deformation potential approximation, the electron- 
phonon coupling is given by ||25|: 



Hei-p = 2 X! J '^^ il^lipsdxd , (2.2) 

where g is the electron-phonon coupling constant (the 
numerical factor •\/7r/2 is introduced to simplify later al- 
gebra). As mentioned in the introduction, alternative 
choices for the phonon interaction are possible, as is the 
case when one studies conducting polymers [^,0 or 
molecular crystals However, for mesoscopic systems 
such as GaAs heterostructures in the metallic regime and 
at low temperatures (T < IK), the acoustic phonons 
which are coupled by the deformation potential to the 
electronic degrees of freedom constitute a realistic model 
of the lowest energy modes of the system. 

We now turn to the interacting electron system. It has 
been shown |p7| that the low energy and long wave length 
behavior of a system of spinless fermions with short range 
interaction can be described by a continuum Hamiltonian 
expressed in terms of bosonic fields. A similar treatment 
exists for fermions with spin We restrict ourselves 

to correlated electron models with nearest neighbor hop- 
ping and short-ranged repulsion/attraction. The field 
operator describing electrons with spin s — ±1/2 (=|, J,) 
is decomposed into a right and left moving part: 

M^) = e^'""V's+(a;) + e-''=-^V.-(x) , (2.3) 

The right (left) moving field operator ips+ (i's-) is then 
expressed as an exponential of bosonic fields: 

ips±{x) = exp{iy^{±ips{x) - es{x))) , (2.4) 

where 9six) = dx'Ils{x')), and the bosonic field Hs 
is the canonical conjugate of ifg. We then introduce the 
bosonic fields describing the charge and spin fields: 

fp^ ^{fl +^1) , (2.5a) 

f'y = ^{'Pl ~ fl) , (2.5b) 
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and their respective canonical conjugates Hp and H^- 
With these definitions the electron density operator be- 
comes II2II: 



dx^p 



(2.6) 



where we have kept only the slow spatial variations of the 
densit y op erator . Ex pressed in terms of the variables in 
Eqs. ( ^.5a ) and ( 2.5b| ), and their canonical conjugate, the 
Hamiltonian describing the interacting electron separates 
into a charge and spin contribution. Each contribution is 
of the sine-Gordon type. Schulz |2^,0 has argued that 
these Hamiltonians can be replaced by "free" , quadratic 
Hamiltonians, provided that the parameters of each con- 
tributions are properly renormalized by the interactions. 
The quadratic Hamiltonian of Ref. | p9[ is our starting 
point: 



with 



Hn + Hr, 



ri,. 



(2.7) 



(2.8) 



with V = p, cr. The free electron gas limit corresponds to 
Kp = K(j = 1 and Up = — vp, the Fermi velocity 
of non-interacting electrons. Note that we have omit- 
ted an additional cosine term |2^ in the spin channel: 
this marginall y irr elevant contribution will be discussed 
briefly in Sec. VB . 

For the Hubbard model with on-site repulsion, the jus- 
tification for this "free" Hamiltonian is based on the 
fact that the large U and the small U limit represent the 
same phase of the model, and there are no singular points 
between these two limits. Consequently the interaction 
renormalize to the "fixed point" Ka- = 1. The remain- 
ing parameters -ftTp, Up and of the Hamiltonian in Eq. 
(2.7) can be obtained by solving numerically the equa- 
tions describing the ground ||3l| and excited states [ p2| , 
as will be shown later in Sec. |v4 For other, more gen- 
eral models where an exact solution is not available, the 
parameters Up, Ua-, Kp, K^j can be obtained, for example, 
by finite size diagonalization for small systems. 

The electron-phonon coupling originates from the elec- 
trostatic interaction between the electrons and the atoms 
of the solid. Consequently, the phonons couple only to 
Hp, since it is the total charge density from up and down 
spins that matters. 



III. LAGRANGIAN FORMULATION 



To calculate thermodynamic quantities, we use a La- 
grangian formulation | |33| where the partition function is 
represented as a functional integral over the phonon (ph) 
displacement field and the bosonic fields associated with 
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the charge and spin degrees of freedom. For example, 
if we wish to calculate the thermodynamic average of a 
functional f[(pv, ^v] which depends on the conjugate op- 
erators (pn and n^, with v = p^a, we can express this 
average in terms of a functional integral as follows 

<: fWu.i^u] ■■>= Z-^ j Dip.DU^DdDUd exp{-5£[^„ H,, d, HJ} 

x/[^,,n,], (3.1) 

where Z is the partition function. The symbol : / : de- 
notes the normal ordered form of f with the momentum 
operators brought to the left of their conjugate field oper- 
ators. We shall see below that this type of normal order- 
ing only leads to oscillatory phase factors and does not 
affect the asymptotic behavior (in the thermodynamic 
limit) which is our main interest. These oscillatory terms, 
however, are relevant for finite size effects [Q. The Eu- 
clidean action Se can be decomposed as follows: 

SE[^^,Tl,^,d,Ild] = dr / dx [Lph + Lp + L„ + Lp^ph] ■ 
Jo Jo 

(3.2) 

Here, r is the imaginary time coordinate (/3^^ ~ ksT). 
The corresponding Lagrangian densities for the phonons, 
the charge and spin degrees of freedom, and the coupling 
term are defined by: 



Lph 



1 



2C 
u,,K, 



[d^df - i\lddrd 



(3.3a) 



-n 



2 " 2K, 
Lp-ph = g{dxd){dx<Pp) , 



iliudripi, (3.3b) 
(3.3c) 



with V = p,a. Note that the Lagrangian density describ- 
ing the system is quadr atic in the fields, and the func- 
tional integral of Eq. (3.1) will be performed exactly. 
This is a direct consequence of the bosonization method, 
as the fermion density is linear in dx^p (see Eq. (2.6)). 



In principle, the functional integrals of Eq. (3.1 



plicitly contains summations over the winding numbers 
associated with the zero- modes of the fields |18 . These 



zero-modes, however, shall be ignored here as we are cal- 
culating only bulk quantities, but they should be included 
in the calculation of any finite size effects. Consequently, 
we shall omit in the Fourier sums below the A; = con- 
tribution (this constraint becomes irrelevant in the ther- 
modynamic limit). 

Th e cro ss terms —illddrd and -~ill^dr(pi, in Eqs. (3.3a) 
and ( pb| ) can be eliminated if we redefine the canonical 
conjugates as follows: 

lid —> lid - iCdrd 



The Lagrangian densities for the phonon, charge, and 
spin degrees of freedom then become: 
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Lvr. = Yc^l + ^{d.df + ^{drdf (3.4) 

(3.5) 

Later on, in the calculation of the Green function and 
the correlation functions, we will have to keep in mind 
this shift in the canonical conjugate 11,^. 

As we shall later be concerned with the calculation of 
one and two particle Green function for electrons only, we 
can integrate out the phonon degrees of freedom at this 
point. We must therefore perform the following Gaussian 
integral: 

h[Vp\ = J Dd exp - J dx{^ d[-dl - c^dl]d 



gddl^pp) 



(3.6) 



where dx = dxdr. Denoting by [—9^—0^9^] ^ the inverse 
of the operator act ing on the phonon field on the right 



hand side of Eq. (3.(:), the functional integration over 



the phonon field gives: 

/i[<^p] =ciexp[|^ j d^d^{dl^,){^)[-dl-^dl]-\^,^)]{dl^,){^) 

(3.7) 

where ci = [det(— 9^ — c^9^)]^^/^. We can calculate the 



propagator in Eq. (3/7) by going to a Fourier repre- 
sentation for the charge field thereby assuming periodic 
boundary conditions in space and imaginary time: 

with k = (k,uj), X = {x,t) and k-x = kx-\-ujt, and where 
the fc = term corresponding to the zero-mode is to be 
omitted henceforth. We then obtain the contribution for 
the effective action for the charge degrees of freedom: 

where we have introduced the inverse propagator for the 
charge field Lpp-. 

where b — g^J KpUpjUC, (with ?i's restored) is a constant 
which has the dimension of a velocity squared. The first 



two terms of Eq. (3.10) represent the contribution of the 
uncoupled charge field. The coupling to phonons there- 
fore creates an additional quadratic contribution to the 
action, which contains all retardation effects: a collec- 
tive excitation of the interacting electron system may 
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create a phonon at time t and reabsorb it at a later 
time t' . We note that k^/DpCk) is the effective charge 
density propagator in Fourier space. This retarded in- 
teraction is always attractive, and will be responsible for 
the formation of Cooper pairs. We also note that this 
retarded interaction term enters with a minus sign ir- 
respective of the sign of the coupling constant g. As a 
result, an instability arises at zero- frequency and when Up 
approaches b/c (from above) because Dp then approaches 
zero, signaling an instability in the charge density propa- 
gator, /Dp{k), towards long wave-length density fluc- 
tuations. This instability has been known for a long time 
and has first been recognized by Wentzel and later 
by Bardeen [ pT| , its consequences for the thermodynamic 
(static) properties of a non-interacting system have been 
discussed in Ref. This instability will be of great 

importance in our following discussion of the dynamic 
properties of the Luttingcr liquid, in particular with re- 
spect to its various possible phases characterized by the 
retarded response functions. 

Had we treated the electrons with conventional fermion 
operators the integration of the phonon variables would 
have resulted in an effective interaction which is quartic 
in the fermion operators. Moreover, in problems where 
a single particle is coupled to a phonon bath |Q, the 
"influence functional" which results from the functional 
integration of the bath variables is in general quite com- 
plicated. The particle is then coupled to an effective po- 
tential which contains "memory effects" [ ^5p^ , and one 
has to take recourse to approximation methods to un- 
derstand how the coupling to the phonon bath affects 
the single particle motion. This is not the case here: a 
one-dimensional electron system (with or without inter- 
actions) coupled to phonons can be treated exactly if one 
uses the bosonization description for the charge and spin 
fields. 



IV. SINGLE PARTICLE PROPERTIES 



In this section we discuss the effects of the electron- 
phonon interaction on the single-particle properties. In 
particular, we show that an arbitrarily small coupling to 
phonons immediately drives the system into a non-Fermi 
liquid state, even in the absence of electron interactions. 



A. Green function 



We proceed to calculate the single-particle Green func- 
tion: 



G,(x,t) = - < T^,(a;,T)v4(0,0) > 



(4.1) 



where T is the imaginary time ordering operator. Using 
the decomposition of the fermion operators into right and 
left moving components (Eq. (2.4), the Green function 
becomes: 



9 



G,(x,r) = - J2 e'"'"" <rV'sa(a;,r)vL(0,0) > , 

(4.2) 

where ipsa (x, t) is expressed in terms of bosonic variables 
as in Eq. ( p.4| ). The diagonal contribution a = a' turns 
out to be the dominant one. We thus need to calculate 
the expression 

X 5„(x,t) , (4.3) 



with the charge and spin fields introduced in Eqs. (2.5a) 



and ( p. 5b ), and wi th = 1, and s^- = s. The last line on 



the r.h.s. of Eq. (O) is a c-number which arises from 
the normal ordering procedure, and, moreover, contains 
the information about the time ordering of the Fermi 
operators: 

5n(::c,T) = + J] e"^[''''(°^o)'^-(^'^)l/2, for r>0 

gn{x,T) = - n e"^!''''^-'^)'^-^"'^)!/^, for r < . 

(4.4) 



The normal ordering contribution, Eq. (4.3), involves the 
commutators between Oi, and tpi, evaluated at unequal 
times. The calculation of these commutators is given 
in Appendix A, where also the explicit time dependence 
of the field operators is derived. The normal ordering 
contribution then becomes: 



gn{x,T) = sgn(T) Jl 



e 

I/— p,cr 

ia-TT 2:sgn(r) / L 



-a7rsgn(r)[e„(3:,r),ip„(0,0)]/2 



= sgn(r)e^"'^^''S"^^^/^ , x ^ Q , (4.5) 

wher e the exp licit form of the commutator is given in 
Eq. ( A8a , A8b| ). For later use we also give the limiting 



result for T — > 0''=, 

where z{x) = exp(27ria;/L — e) (eL is a short- distance 
cutoff of the order of the lattice constant). Eq. (|4.6|) is 
valid for arbitrary x. 



Next we compute the expectation value in Eq. (1.3). 
Going to the Fourier representation, 

n-W = 7TrE^^*'^''n.(k): (4-7) 

PL ^ 

we perform the integral over the canonical momenta Hp, 
lia with the result: 
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< exp 



exp 



> 

ik-x _ 2|2 



fc2 



(4.8) 



Note that the exponent on the r.h.s. ofEq. ( |4.8| ) contams 
a negative infinite constant due to the summation over 
Matsubara frequencies, however, as we shall see below, 
this contribution is exactly cancelled when integrating 
out the field tp,^. 

The functional integral over ipi, gives the contribution: 



/i/(x) =< exp 



2 

TT 1 



2 u^Ki, 





dx' dr(pv{x' , t) 



dx' dr'~Pv{x' , 0) 



> 



(4.9) 



where the last term in the exponential comes from the 



shift in the conjugate momentum as discussed in Sec. Ill 
The propagator for the spin field tp^ is given by 



Using this result and Eq. ( [3.1C| ) we obtain for ly 



(4.10) 



/j,(x) = exp 



4/3L 



k 

,zk-x 



- 1 



4iasi/u; 

kUjyKy 



cos(k ■ x) 



(4.11) 



The first term in the exponent of Eq. (4.11) diverges 
due to the sum over Matsubara frequencies. This term 
is positive, and the divergent part is cancelled by the 
contribution of Eq. (4.8). To see this it is useful to 
write: 



to 



fc2Z),(k) 



(4.12) 



When inserted in Eq. (4.11), the second term on the 
r ight hand side cancels exactly the contribution of Eq. 

(S- . . ^ 

The thermodynamic average in Eq. (4.3) then be- 
comes: 

<T^.„(x)VL(0)>-;^ n -i^E^-~'W 
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u,,K,,k 



sin(fea;) sin(a;T^ 



(4.13) 



Note that the exponent in Eq. ( 4.13| ) has a real and imag- 
inary part coming from both the thermodynamic average 
and the normal ordering procedure. This implies that the 
Green function has an oscillatory component and decays 
with a power law as e xpect ed. The frequency and mo- 
mentum sums in Eq. ( 4.13 ) are performed in Appendix 
B. The result is: 



<TV'.„(x)V'L(0) >=^exp 



1 .Cg 

4 '^^Kr, 

f3=± P 



^^Fi3)\n[\l - z{x + ivfi\i 



-\{K, + ^)H\l^zix + tu,\T\)\/e] 
sgn(r) In ~ *w<t|t|) 

4 1 — z{x + iUcr\T\) 

-sgn(r)— > Frtln- -. . . ., 

" l-zix + tVfs\T\) 

X g„(x) . (4.14a) 

We have introduced effective velocities induced by the 
phonon coupling: 



as well as the abbreviations 



Up vl-c^ 



v± 

o,2 



(4.15) 

(4.16a) 
(4.16b) 



The thermodynamic limit corresponds to taking \x + 
ivT\/L << 1, which in turn implies l—z{y) ~i2Try / L+ 
e. We can now give the final expression for the Green 
function in t he thermodynamic limit using Eq. ( pla| ) 
and Eq. (EJ): 



G«(x,r) 



1 

eL 



eL/2TT 



xn 

X 

Q = ± 

7=± 



X + lUaT 

eL/2TT 



KpUpF^/iVfi+Cp/iKp 

X + ivpT 

, ■ I I \ sqn(T)/2 

sax + iU(j|t| \ ^ ^ " 



l\x + lUaT 

ax -I- iv^\t\ ^ ' ^' 
i\x + iv^T\ 



9n{x,T) 



(4.17) 



Near x = one has to add eL/2-K to v\t\. At large 
distances and/or times (?„ becomes unity and the power 
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law decay of the Green function becomes then (in obvious 
notation) 



Gs{^)^\x,t\-^-' , (4.18) 

with 

4 4Ko- 4Ap 4 ^ ^ 

and the electron-phonon parameters A and B are defined 

by 



/3=± 

(4.20a) 

B{up, c,b)^J2^0 = fl + ^) (4-20b) 

These exponents play an important role in the following 
discussion of the ordering fluctuations. For later use we 
only note here their limiting behavior. Defining a critical 
charge velocity by the equality b/upC — 1, or alterna- 
tively: 

- (4.21) 



k; ncc 



we see from Eq. (4.15) that tends to zero, and VZ- to 



u*p^ + c^, as Up approaches this critical velocity u* from 
above. As a consequence, the exponent A diverges to in- 



finity and B decreases to the finite value u*p/ \J u*"^ + (? as 

Up u*. For Up < u* the velocities become complex and 
the model becomes unphysical. Thus we must require 
that Up/Kp > u*/K*, or equivalently that b/{cup) < 1, 
the equality sign defines the Wentzel-Bardeen singular- 



ity |37|. We emphasize that this singular behavior is a 
non-perturbative effect of the electron-phonon coupling 
and originates from the instability of the effective elec- 



tron propagator as pointed out after Eq. ( [3.10 ) . We also 



note here that the limits 5 — s- and Up — > u* do not 
commute. 

In the limit of vanishing phonon coupling constant, the 
charge and spin propagator have an identical structure, 
and A, B 1. The result is: 

hm G,(x) oc la;, r\-^./i-UiK,-Kj4-i/4K, (4 22) 

In particular, for non-interacting electrons Kp — — 
1, we recover the 1/x power law dependence of the Green 
function of the free electron gas. 

To illustrate our results with a specific model of inter- 
acting electrons, we plot the quantity 5 for the Hubbard 
model in Fig. |l|a and |^b. The Luttinger liquid parame- 
ters Kp, Up are determined numerically later on in Sec. 



VI 
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We first plot i5 for a quarter filled band as a function 
of U in Fig. |l|a, for several values of the phonon coupling 
parameter. For small phonon coupling, S increases mono- 
tonically with increasing U. As the phonon coupling is 
further increased, S acquires a minimum. 

Next, we choose U = 2 and vary the filling factor from 
zero to half filling in Fig. |^b. 6 increases dramatically 
when these two extremes are reached, as the correlation 
effects between electrons dominate the physics in both 
cases. In particular, (5 — > oo at the Wentzel-Bardeen 
singularity. 



B. Momentum distribution function 



We now examine the effect of the phonon coupling on 
the momentum distribution function which is given by 



1 /■+°° 

^ „ J-oo 



(4.23) 



Next, we insert the result of Eqs. ( |4.18| ) and ( 4.19 ) at r = 
0^ Pq] , use that in the thermody nam ic limit gn{x, 0^) = 



ibiala; + ia\/{x + ia) (see Eq. ([l.q)), and find (setting 
Oo = ei/27r), 

^ 1 ^ ° • (4.24) 

2m \x -\- iao\ X A- laao 

As pointed out in Ref. p9[, for a system with (5 > 0, the 



exponential in the integrand of Eq. (4.24) is not neces- 
sary any more to insure the convergence of the integral. 
For k ~ akp^ we may omit this oscillatory factor, and 
obtain the distribution function evaluated aX k = akp 

- r(i/2 + V2) (4,3) 



2V^r(l + ,5/2) ' 

where r(a;) denotes the Gamma function. In the limit 
(5^0, Na{akp) = 1/2, which differs from the Fermi 
step function obtained for the free system at 5 = 0. The 
momentum distribution function has a finite value at the 
Fermi wave vector, but its derivatives are all singular at 
that point (see below). This result was previously found 
in Ref. | |39| ] for a system of interacting fermions but with- 
out phonons. The fact that the momentum distribution 
function has no jump at the Fermi surface for interacting 
fermions expresses the failure of the Fermi liquid quasi- 
particle picture to capture the low temperature behavior 
of interacting one-dimensional systems, which has been 
known for the Luttinger model since the pioneering work 
of Mattis and Lieb Q . 

We can further determine the behavior of the momen- 
tum distribution function around k = kp hy noting that 
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N+{k) = - 



dx 



cos {aox{k — kp)) 

(x2 + l)l+'5/2 



■{l+0{^\k-kF\)) 



Evaluating the tabulated integral we then find 
N{k ~ kp) ^ N+ikp) - csgn(fc - kF)\k - kp]^ 



(4.26) 



(4.27) 



with c^ai ^/^r(-l - (5)/r(-(5/2). Eq. (|^ exhibits 
the singular dependence of the momentum distribution 
function at k — kp. Note that for the pure electron 
system (g ^ 0), taking = 1, we recover the exponent 
of the momentum distribution function given in Ref. [ p9| , 
namely: 



1 



(4.28) 



C. Free electron gas coupled to phonons 



In Eq. ( 4.19 ), both the effect of the interaction between 
electrons and the effect of the coupling with the phonons 
are included. In this section, we specialize to the case 
of free electrons coupled to phonons in order to clarify 
the role played by the phonons. In this case, S takes the 
form: 

S ^ A/4: + B/i- 1/2 . (4.29) 
It is convenient to rewrite A in the following way: 

A = B + f ^ , . (4.30) 

Noting that 

= l-by[up{v++v-)]^ <1, (4.31) 

we calculate 

b \V&2 



\v+v-(v++v-)j \uf, 



(4.32) 



With Eq. ( py2q ), and A, B > this imphes that 6 > 0, 
for g ^ 0. If we start from a free Fermi gas, and couple 
it to phonons, the system is driven away from the free 
Fermi behavior and the Fermi distribution function is 
totally destroyed for an arbitrarily small phonon coupling 
strength g. 



V. ORDERING FLUCTUATIONS IN 
ONE-DIMENSION 

In this section, we determine which type of fluctua- 
tions dominate in an interacting electron gas coupled to 
phonons. 
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A. Correlation functions 



The tendency of the system towards ordering mani- 
fests itself in divergent correlation functions, which de- 
scribe the response of the system to external perturba- 
tions. Thus, to characterize the type of long-range order 
(collective state), which might be realized in higher di- 
mensions, we follow standard practice [3P,p9[| and con- 
sider response functions of the type 



R{x,t)=-ie{t)< [0''{x,t),0] >, 



(5.1) 



where the charge-density wave (CDW) response is gen- 
erated by the operator, Ocdw = [V'+t^'-T + V'+x'/'-i]/2, 
the transverse spin-density wave (SDW) response by 
OsDW = tp+s'^--s, singlet pairing (SS) by Oss = 
tl>+s'^--s, and triplet pairing (TS) by Ots = tp+si^-^s- 
The corresponding two-particle Green functions in Mat- 
subara representation read then: 

N{x,t) = -ie^^feFi; < TV^l^(a;,T)V'+T(a;,T)V'||V^-T > 

(5.2a) 

-ze'2fc-- < T^1^(x,t)V'+t(x,t)^;^^_X > 

(5.2b) 

-i < TV'-i(a;,T)V'+T(2:,T)V^.|^/'lx > (5.2c) 
-i < TV'_-r(x,T)V'+T(a;,T)V'+|i/'l| > , (5. 2d) 



As{x,t) = 



with the notation -0^3 defined in Eq. (2.4). N repre- 
sents the correlation function associated with the CDW 
state, which for example occurs when the nearest neigh- 
bor repulsion is large compared to the other parameters 
of the model, x represents the transverse SDW response 
function, which describes the onset of antiferromagnetic 
ordering. As (At) gives the probability amplitude for 
a singlet (triplet) Cooper pair emitted at the origin of 
space and time to be annihilated at x and r, and thus 
characterizes the superconducting fluctuations. To sim- 
plify notation we have retained only a representative term 
in the correlation functions, and this term will determine 
the correct small frequency /momentum behavior of the 
whole correlation function. 

Using the bosonized version of V'±s in conjunction with 
t he de composition into charge and spin of Eqs. ( 2.5a ) and 
( ^.5b| ), we obtain (for r > 0): 



i2kf 



N{x,t) = -'- 



{eLf 



< 



n 



exp 



«V2^(^,(a;,r)-(^,(0,0)) 



> 



As(x,r) 



{eLf 



{eLf 



< exp 



exp 



< exp 



(5.3a) 



-iV2^i dx'n^{x',T)- dx'll„{x\0) 



lV2TTi^pix,T)~ipp{0,0)) 



> 



(5.3b) 



-,-v^(/..'n,(.',.)-/..'n,(.',o)^ 
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exp 



At(a;,r) 



(eL) 



2 < exp 



> 



(5 



-iV2^\^l dx'np{x',T)- I dx'Up{x',0) 



exp 



-iV2^(^J dx'n^{x',T) J dx'U^{x', 



0) 



Th e fu nction al integrations which are imphed in Eqs 
(5.3a) to (5.3d) have already been encountered in the cal- 



culation of the single-particle Green function in Sec. IV 



We therefore omit the details of this calculation. Also, 
because we are only interested in the long distance and 
long time properties of the correlation functions, we have 
ignored the oscillatory terms arising from the normal or- 
dering procedure. The results are: 



N{x,t) 



le 



i2kj='X 



{eLf 



■ exp 



Kp J2 —Fp In |1 - z{x + ivpT)\/e 



-K„ In |1 — z{x + iVaT)\/e 



(5.4a) 



■ exp 



'KpY^ ^F0\n\l - z{x + ivpT)\/e 



--— In |1 - z{x + iUaT)\/e 



(5.4b) 



As(a;,r) 



■ exp 



/3=± 

Kc In |1 — z{x + mo-T)|/e 



(5.4c) 



At(a;,r) 



■ exp 



--— In |1 — z{x + iUaT)\/e 



(5.4d) 



with C± and F± defined in Eqs. ( [4.16a|) and ( [f.l6b|) . 

The thermodynamic limit is then taken to obtain the 
power law behavior of the correlation functions: 

N{x, t) cx e'^'^^^ n 1^ + ivpT\-''^^''^^'''(' \x + iu^T\-^'' 
/3=± 

(5.5a) 

X(x, r) (X e'^'^^^ W \x + ivpT\-''''^''^i'/''('\x + iUaT\-^''^'' 
f}=± 

(5.5b) 

As(x,r)(x \x + ivfjT\^^'i'^'^''\x + iu^t]^^" (5.5c) 
P=± 

At{x,T)(x Y[ \x + ivjiT\-^i'l^''\x^'iu„T\-^l^' . (5.5d) 
/3=± 



(5 
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B. Divergence criteria 



The indication that a given type of flucuations domi- 
nates in the system is provided by the divergence of the 
Fourier transform of the corresponding correlation func- 
tion at low frequencies and small momentum relative to 
q = 2kF {q = 0) for N an d x (A^ and At). 

Inspecting Eqs. (5.5a) to ( ^.5d ), we notice that such 



a divergence can only originate from the power law be- 
havior of the correlation functions at large distances and 
imaginary times. Whether or not there is a divergence 
then depends on the exponents which were calculated in 
the previous section. In particular, we note that 

\X + iVmaxT\~'"' < \Ck{x,T)\ < -|- iu„i„T| (5.6a) 



whereCfe is given by the corresponding rhs of Eqs. 5a ) 



to ( |5.5d[ ), 7fe by the associated sum of exponents, and 
Vmax/min = mux / min{v+, V- , Ua} ■ From this, we ob- 
tain immediately the divergence criteria for the Fourier 
transforms of the correlation functions. The condition 
for CDW fluctuations to occur is then, from Eq. ( |3.5a[ ): 

AKp + K„ < 2 . (CDW) (5.7) 



From Eq. ( 5.5b ), SDW fluctuations occur when 
1 



AKp + — < 2 . (SDW) (5. 



From Eq. (5.5c), singlet Cooper pairs fluctuations are 



present in the system when: 

-^+K,<2. (singlet) (5.9) 

Finally, we will have triplet superconducting fluctuations 
when: 

13 1 

— + — < 2 . (triplet) (5.10) 

Ap Act 

For models with interactions which do not depend on 
the spin {SU{2) symmetry), such as the Hubbard model, 
K„ = 1, and the criteria for CDW and for SDW or- 
der (for singlet and for triplet pairing) are a priori the 
same. However we expect SDW and triplet Cooper pairs 
to dominate: CDW order will not occur because SDW or- 
der builds up before the longer range type of fluctuations 
have time to settle in the system Q. Similarly, triplet 
superconductivity will dominate over singlet supercon- 
duct ivit y. Th is is not noticeable in the criteria of Eqs. 
( ^.7[ ) to ( ^.10 ) because additional logarithmic corrections 



have been neglected here For the electronic models 
we are considering in this paper, these corrections arise 
from marginally irrelevant operators in the spin chan- 
nel, which effectively generate longer tails for the SDW 
(triplet) correlation function than for the CDW (singlet) 
correlat ion functions. These operators have been omitted 
in Eq. (|2.8|) for simplicity. If both divergence criteria for 
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spin waves and triplet Cooper pairs are met for a given set 
of parameters, the correlation function "which diverges 
the most" determin es th e phas e of the system. Note that 



the criteria of Eqs. (5.7)-( p.l0 ) are valid for any interact 



ing electron system with a finite range interaction: what 
remains to be done is to study the dependence of Kp and 
Up on the interaction parameter(s) of the electron Hamil- 
tonian and the filling factor for the electrons. Note that 
for repulsive interactions, we have Kp < 1. 

At Kp = Ka = 1, the criteria for SDW order and 
triplet pairing take the simple form: 

A<1 . (SDW) (5.11a) 
B < 1 . (Cooper pairs) (5.11b) 



from the results of Sec. IV C and noting that 



1 — b'^ /(cUpY 

with Up < v+ + V- and c > 0, we conclude that ^ > 1, 



and Eq. (5.11a) is never satisfied , wh ile Eq. (5.11b) is 



always satisfied according to Eq. ( 4.31 ). The Cooper in- 
stability is always present for a system of non-interacting 
electrons, regardless of the magnitude and sign of the 
phonon coupling constant. 

When no phonons are present in the system, A = B = 
1, and a system of electrons with repulsive on-site in- 
teraction alw ays has Kp < 1, which implies from the in- 
equality (U) that SDW (CDW) fiuctuations dominate 
over superconducting fluctuations, as expected. 

A check of the above results can be obtained by consid- 
ering non-interacting electrons (Kp = K^ = 1) coupled 
to phonons in the limit where the speed of sound is much 
larger than the Fermi velocity. If we set c — > oo, keep- 
ing^/ cUp finite (0 < hjcUp < 1), the propagator of Eq. 
( |3.10 ) becomes: 



Dp{k) ~ {KpUp)-\Lj^ + u^e - h^e/c^ + c-^Oiiv^k'^)) . 

(5.13) 

This propagator now corresponds to a system of elec- 
trons (with no phonons) with an attractive interaction 
without retardation effects, as in the original Cooper 
problem. We should therefore be able to reproduce the 
result of Luther and Peschel for spinless fermions 
which predicts strong pairing fluctuations for the attrac- 
tive case. In the limit c — > oo, we obtain from Eqs. 
( pOal ) and the result that A = (1- (fe/cMp)^)"^/^ 



and B = {1 — {b/cUpYY^^ . The pairing correlation func- 
tion for triplet pairing diverges, according to Eq. ( ^.10 ) 



but the criterion for SDW order is not met. The at- 
tractive interaction induces strong pairing fluctuations, 
in anology with the Cooper problem and the result of 
Luther and Peschel for spinless fermions. 
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VI. HUBBARD MODEL 



In this section, we shall consider the effect of phonons 
on the correlation function exponents derived in the pre- 
vious sections, for specific interacting electron models. 
Most of the discussion will consider a Hubbard model 
with on-site repulsion parameter U: 

Hel = -tj^i'is^'^^+i,^ + + U ^71,1^,1 , (6.1) 

i.s i 

with t the hopping matrix element. We shall also con- 
sider the extended Hubbard model with nearest neighbor 
repulsion parameter V in the limit [/ — )■ cxd, where ana- 
lytical expressions for Up and Kp can be found at quarter 
filling. 



A. Weak interactions 



We first consider the on-site parameter as a pertur- 
bation. The charge velocity is then given by the Fermi 
velocity: 



vp — 2tsm{mr/2) 



(6.2) 



for a system which has n electrons per site (half filling 
corresponds to n = 1). The shift in energy (per site) due 
to the on-site repulsion term in Eq.(6J) is then AE = 
Un? /A:. One therefore deduces the first order expression 
for K„: 



1 - U/ttvf + OiU"^) 



(6.3) 



Inserting the latter expression in the criteria for SDW 
order and triplet pairing, we obtain: 



2U 9{c^ + 2v+v-) 



ttvf 
2U 



62 



(SDW) 

(Cooper pairs) 



(6.4a) 
(6.4b) 



with v± given by Eq. (4.15). Eqs. (6.4a) and (6.41:) dis- 
plays the effects mentioned in the Introduction: the cou- 
pling to phonons will induce Cooper pairs if the coupling 
constant is large enough to overcome the instantaneous 
repulsion between the electrons. Similarly, in order to 
have SDW fluctuations, the on-site repulsion parameter 
must be large enough to overcome the retarded electron 
atttractive interaction mediated by the phonons. For suf- 
ficiently small phonon coupling v+ —> c and v- —^ vf, 
the above inequalities are mutually excusive: there is a 
region in parameter space: 



VpivF 



2U P{1 + 2vf/c) 

TTVF Vp{vF + c)2 



(6.5) 



which separates the SDW order and the triplet pairing 
phase. 
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Note that the perturbation expansion breaks down in 
both hmits n — > an d n — > 1. In these two cases, we 
shall see in Sec. VI B that Kp — > 0.5, which signals that 
the correlation effects are important. Near n = 1, Umk- 
lapp scattering becomes relevant, and the system feels 
the proximity of the metal-insulator transition. 

We plot the phase diagram of the system for [/ — 0.01 
in Fig. H, using the perturbative result of Eq. (6J), as a 
function of the filling factor n and the phonon coupling 
parameter b/cup. We choose c eLt/h, which is reason- 
able for heterostructure systems where c ~ 5 x 10^[cm/s] 
and vp ~ 10^. The superconducting phase is separated 
from the SDW phase by a metallic regi on, as suggested 
by the small coupling result of Eq. (|6.5D. 



B. Hubbard model at arbitrary filling factor 

To determine the Luttinger liquid parameters Kp and 
Up for arbitrary U and arbitrary filling factor n, we fol- 
low the method proposed by Ref . . We calculate the 
ground state energy by solving numerically the integral 
equation of Licb and Wu [§1^,0 : 



f{k) = ^ + -cos(fc) j'^ dfc'i?^-(sin(/c) -sin(fc'))^,/(fc') 

(6.6) 



27r u 



where u = U /t, and f{k) is a distribution function eval- 
uated at pseudo-momentum k. k take values in the in- 
terval [— Q,(5], and Q < TT depends on the filling factor 
through the relation dkf{k) = n. The function R(x) 
is the cosine Fourier transform of 1/(1 -I- e^) evaluated 
at x/2. With these notations, the ground state energy is 
given by: 

rQ 

EG = -2t dkcos{k)f{k) . (6.7) 

A first relation between Kp and Up is obtained by noting 
i9|,i0|: 

d^Ea _ nup 

The velocity of the charge excitations is obtained from 
the particle (or hole) excitation spectrum, which was cal- 
culated by Coll . Up corresponds to the group velocity 
of particles (holes) in the limit of small momentum: 



Ae(fc) 



Q rQ 



Up = lim — )tt = 2sin(Q) - fi / dkg(k) - 2t / dkcos(k)g(k) /27r/(Q) 
fc^o Ap{k) \ J_Q J 

(6.9) 



with Ae(fc) the particle (hole) excitation energy, Ap{k) 
its associated momentum a nd u the chemical potential. 
The function g{k) in Eq. ( |6.9| ) is the derivative of the 
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particle (hole) excitation distribution function in pseudo 
momentum space. It is the solution of the integral equa- 
tion: 

g(k) = ^ cos(fc) cos(Q)_R' I — (sin(A:) — sinlQ)) 

+ ^ cos(fc) p dk'R(^^{sm{k) - sin(fc'))^ g{k') , 

(6.10) 

where R'{x) denotes the derivative of R{x). 

Once the parameters Up and Kp are specified, we can 
deter mine i n wh at region of parameter space the criterion 



( p. 10]) and ( p.SD are satisfied. The spee d of sound is as- 



sumed to have the same value as in Sec. VIA. We choose 
to represent the phonon coupling with the dimensionless 
parameter b/cup. We emphasize that: 



— (6-11) 

CUp c y TiQup 

is an effective electron-phonon parameter, which, strictly 
speaking, depends of the filling factor. Nevertheless, it 
allows us to display the results in an obvious fashion. 
First we consider the case of h/cup strictly smaller than 
unity (i.e. excluding the Wentzel-Bardeen singularity). 

First, we vary the filling factor n and the phonon cou- 
pling for fixed U . Our results are illustrated in Figs. 
(|a)-(|c), for U/t = 0.7, U/t = 0.3, U/t = 0.1 . The re- 
gions where SDW fluctuations dominate, or alternatively, 
superconducting fluctuations dominate, are se parat ed by 



a metallic region, as was proposed in Sec. VIA . For 
small phonon coupling, SDW fluctuations dominate for 
arbitrary filling factor. As the coupling is increased, we 
cross the metallic region first near quarter filling, and 
eventually we reach the regime of superconducting fluc- 
tuations. For low filling factors and near half filling, 
correlation effects dominate the picture, and supercon- 
ducting fluctuations cannot be sustained, even at large 
phonon coupling. By comparing Fig. (||a) to Fig. (||b) 
and Fig. (||c), we observe that the region of parame- 
ter space where superconducting fluctuations dominate 
shrinks as U is increased. This demonstrates the inter- 
play between the instantaneous repulsion and the phonon 
mediated attraction for determining the phase of the sys- 
tem. We warn the reader that due to some limitations in 
numerical accuracy associated with the sharp variations 
of the charge velocity in the immediate neighborhood of 
n = 1, Figs, ^a to seem to indicate that the SDW 
phase dominates in this limit. This is not the case, for as 



we shall see in Sec. VIC, antifcrromagnctic fluctuations 
are always suppressed for h/ cUp — 1. 

Next, we choose the filling factor at quarter filling 
n ~ 1/2), and vary U and the phonon coupling, see Fig. 

Again, the metallic region separates the superconduct- 
ing phase from the SDW phase. At low U {U/t < 0.1), 
Kp « 1, and the Cooper instability pushes the system in 
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the superconducting phase. As U is further increased, the 
phase boundaries appear, and the phonon-mediated at- 
tractive interaction is quickly overcome by the instanta- 
neous repulsion between electrons. Beyond U/t ^ 2, the 
SDW fluctuations dominate for all values of the phonon 
coupling. 



C. The Wentzel— Bardeen singularity 



In the preceding section, we plotted the phase diagram 
of the system as a function of the filling factor and the 
on site repulsion versus the phonon coupling parameter 
b/upC. The maximal value for this latter parameter in 
Figs. H and ^ is 1, and corresponds to the Wentzel- 
Bardeen singularity where the compressibility of the sys- 
tem becomes negative. We now show that for the Hub- 
bard model, this singular point is accessible by varying 
the filling factor and discuss the physics in this region of 
parameter space. 

We first reexpress the phonon coupling parameter in 
terms of the bare phonon coupling constant, the charge 
veloci ty, and the parameter Kp, as is apparent in Eq. 
( 3.11 ). Next, we plot the ratio Up/Kp in Fig. || for sev- 
eral values of J7, as a function of the filling factor. Up/Kp 
vanishes both for zero filling and half filling, and in the 
latter case, the variation of UpjKp is most dramatic for 
low V . If we now consider the phonon coupling constant 
to be fixed at an arbitrary non-zero value, we can get 
arbitrarily close to h/upC = 1 by, for example, increasing 
the fi lling factor towards n — 1. As pointed out in Sec. 
IV A, when Up/Kp approac hes f rom above the critical 
value u*p/K* defined in Eq. ( 4.21 ), the quantity A which 
appears in the Green function exponent and in the SDW 
and CDW correlation function exponents diverges^ and B 
decreases to a finite value B* = u*p/ (m*^ -I- c^)"'^/^. Conse- 
quently, SDW (CDW) fluctuations are suppressed as we 
approach the Wentzel-Bardeen singularity, and the sys- 
tem is pushed towards the triplet superconducting phase, 
via the metallic phase depicted in Figs. |^ and ^ Given 
the fact that B* cx and 1/2 < Kp < 1 the supercon- 
ductivity criterion, B < Kp, can be met for sufficiently 
small g ^ 0, by approaching the Wentzel-Bardeen singu- 
larity. The system can be driven in the superconducting 
phase by approaching half filling. 



D. Extended Hubbard model in the limit U ^ oo 

We now add a nearest neighbor interaction term: 

l/^n,n,+i, (6.12) 

i 

to the electronic Hamiltonian of Eq. ( |6.lD and let U 
00. Each site can be occupied at most by one electron, 
so that the system reduces effectively to spinless fermions 
with the substitution kp ^ 2kp. At half filling (which 
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corresponds to n = 1/2 here as the Fermi momentum 
is 2k p), the spinless fermion Hamihonian with nearest 
neighbor interaction can be mapped to the Hamihonian 
of the XXZ model for a spin chain. The latter Hamil- 
tonian can be solved exactly with the Bethe ansatz . 
For this reason, it is possible [{30|j42| ] to get analytic ex- 
pressions for the parameters Kp and Up: 

= (2 + 4sin"i(u)/7r)-i (6.13a) 
Up = TTt\/ 1 — jcos^^v , (6.13b) 

for V = V/2\t\ < 1. The case v > 1 represents a dimer- 
ized insulating phase One can therefore study the 

behavior of the correlation function exponents as the in- 
sulating phase is approached. 

The phase diagram is plotted in Fig. ^. Obviously, the 
U oo limit pushes the system towards the SDW phase, 
but at relatively low V/2t, a strong coupling to phonons 
b/cup « 0.94 destroys the SDW fluctuations. As V/2t 
increases from to 1, the effects of the phonons become 
less pronounced, and the metallic region gradually disap- 
pears. 

VII. APPLICATION TO ELECTRON WAVE 
GUIDES 

The results described in the previous section can be 
applied in a different context. The contribution of the 
Hamiltonian in Sec. ^ which describes the coupling 
between the charge and phonon fields can alternatively 
represent the coupling between two charge fields associ- 
ated with two interacting quasi one-dimensional electron 
strings corresponding to two transverse modes in an elec- 
tron wave guide. 

The question of the stability of coupled one- 
dimensional interacting electrons systems has been the 
focus of recent interest |3|,||,||. Ref. @ argues that 
the coupled system is stable, in the sense that it stays in 
the Luttinger liquid phase, as long as the hopping matrix 
element between the two chains remains below a critical 
value. Alternatively, a renormalization group argument 
proposed in Ref. |^ seems to imply that the coupled sys- 
tem is unstable for an arbitrarily small coupling. Both 
Ref. Q and Ref. Q assume a transverse hopping ma- 
trix element as the mechanism for destabilisation. 

We reexamine this question from a slightly different 
starting point, considering the electrostatic effects be- 
tween the two chains. We will show below that a thresh- 
old coupling will induce superconducting fluctuations in 
one of the chains, and in both chains as the coupling is 
further increased. After all, if the two electron chains 
can be associated with two orthogonal transverse states, 
it is reasonable to assume that in the absence of impuri- 
ties, the (screened) Coulomb interaction between the two 
chains is the basic interaction mechanism which breaks 
the orthogonality between the two transverse states. In 



24 



this situation, Eq. (2.2) represents a local density- 
density interaction: 



Xx)lljis{x)'ljjl^,{x)ll;2s'ix) 



g- I dx{drippi){d^(pp2) , 



(7.1) 



where ipis is the electron field operator associated with 
an electron on chain i with spin s, and tppi is the charge 
bosonic fiel d of the same chain. Note that the last equal- 
ity of Eq. (7.1) follows from the fact that the two trans- 
verse modes have distinct Fermi velocities. Hence the 
fast varying components of the density operators in the 
first line of Eq. (7.1) do not compensate each other as 
in the coupling between electrons and 2kp phonons. In 
the case of quasi-one-dimensional systems created in a 
GaAs/AlGaAs heterostructure, the assumption of a short 
range interaction can be justified because the surround- 
ing metallic gates screen all interactions. The interaction 
potential of Eq. (7.1) is precisely the same form of cou- 
pling which was recently used to study the transport 
properties of arrays of one-dimensional chains. This cou- 
pling describes two effects: first, the fact that the local 
"site" energy of, say, the first wire (1) is raised (or low- 
ered) in the presence of a nearby density perturbation in 
the second wire (2). Secondly, Eq. (7.1) allows electrons 
from 1 and 2 to exchange as a pair. 

The remaining part of the Luttinger liquid Hamilto- 
nian which describes the uncoupled system is then: 



Hel — Hpi + Hcl + Hp2 + Hc2 

with, for v — p,a and i = 1, 2 



(7.2) 



Hi,, 



dx 



(7.3) 



with n,yi the canonical conjugate of ip^i. In general, the 
charge velocities associated with the two electron chains 
can be different. 



Using the results of Sec. Ill, the effective propagator 
describing the evolution in a given chain i becomes: 



£)pi(k) = (KpiUpi) 



(7.4) 



with z ^ j and — g^KpiK p2UpiUp2'h^ ■ 

The Green function exponent associated with electrons 
from chain i then reads: 



4 



1 



4i^„ 



AK 



in 



. (7.5) 



where 
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1 , , o , . 1 



v±^^ (w^i +u%)± -^{u% ~u%f + 462 . (7.6c) 

The coupling between the two chains will favor super- 
conducting fluctuations in each chain, in analogy with 
the phonon case. The correlations functions which pro- 
vide the signature for these fluctuations are read directly 
from Eqs. (5.2c) and ( p3.2d ), with all fermions operators 



belonging to the same chain. 

The criterion indicating strong superconducting fluc- 
tuations i n a g iven chain i is then translated from Eqs. 



and (|5.10D 



< 2 . (singlet) (7.7a) 

B, 1 



Kp, 



< 2 . (triplet) (7.7b) 



It is important to note that this effect does not depend 
on whether the interaction between the electron chains 
is positive or negative. As a result, the strong super- 
conducting fluctuations induced by the interchain cou- 
pling could be observed in wires where only two channels 
propagate, or alternatively in one-dimensional coupled 
electron-hole systems. We can also determine from our 
analysis of the phonon case that the chain with the lower 
charge velocity will be the first to go superconducting. 

For completeness, it should be mentioned that there 
are two additional correlation functions describing super- 
conducting fluctuations for this particular system. These 
occur when a singlet or triplet Cooper pair is formed with 
one electron from chain 1 and one electron from chain 2. 
We shall ignore this additional super cond uctivi ty me cha- 



nism here, because the inequalities (7.7a) and (7.7b) are 
in general satisfied for a lower threshold phonon coupling. 
For specificity, we consider the case of two Hubbard 



chains coupled by the interaction of Eq. ( [7.11) : we refer 



the reader to Sec. VI, where the dimensionless coupling 
parameter has to be replaced by h/upiUp2. Away from 
zero filling and half filling, we see that for a reasonably 
small U /t, superconducting fluctuations will be induced 
above a threshold coupling. Our result therefore seems 
to be in qualitative agreement with Ref. |Q, where the 
two chains, if initially in the Luttinger liquid regime, will 
remain in this regime below the threshold coupling. Also, 
note that in analogy to the electron-phonon case, the 
Wentzel-Bardeen singularity can be exploited to drive 
the system in the superconducting phase. 



VIII. SUMMARY AND CONCLUSION 

This work has focused on the low temperature prop- 
erties of a system of interacting electrons coupled to 
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phonons in one dimension. The low energy properties 
of interacting electrons can then be described with the 
Luttinger liquid picture introduced by Haldane Q for 
spinless fermions, and later on generalized by Schulz ||29| ] 
for electrons with spin. When such a system is coupled 
to low energy phonons with a deformation potential cou- 
pling, the phonons can be integrated out, allowing for an 
exact calculation of the correlation function exponents as 
a function of the Luttinger liquid parameters, the phonon 
coupling constant, and the sound velocity. 

We first considered single particle properties, calculat- 
ing the Green function "dressed" by the phonon modes. 
For a gas of free electrons, the coupling to phonons alone 
is sufficient to introduce a singularity of the momentum 
distribution function at the Fermi surface, as is the case 
for fermions in the presence of an instantaneous, attrac- 
tive interaction: the power law decay of the Green func- 
tion at large distances has an exponent 1 + 6 which is al- 
ways larger than 1. For increasing, but small phonon cou- 
pling, the repulsive interaction further accentuates this 
effect. However, beyond a threshold phonon coupling, 5 
has a minimum as a function of U . This indicates that 
for strong enough phonon coupling, at low U, the instan- 
taneous repulsion and the retarded attractive interaction 
work in opposite ways. 

Next, we considered the many body aspects of the sys- 
tem. In two and three dimensions it is understood that 
the parameters of the electron Hamiltonian can drive 
the system in a definite "phase" at vanishing temper- 
ature, such as a SDW or CDW phase. In the presence 
of phonons, a natural question is to ask whether the re- 
tarded attraction mediated by the phonons can flip the 
system into a triplet superconducting phase. In 2D and 
3D, the interplay between retarded attractive interac- 
tion and instantaneous repulsive interaction cannot be 
accounted for in a precise matter. After all, the BCS 
theory of superconductivity assumes from the start 
an attractive interaction without retardation effects. The 
interplay between attractive retarded interaction and in- 
stantaneous repulsive interaction can be resolved exactly 
in one dimension. A set of inequalities was derived, which 
determine which type of fluctuations are present in the 
system. 

To illustrate our method, we applied our results to the 
Hubbard model, where the Luttinger liquid parameters 
can be determined through a numerical computation, 
for arbitrary filling factor and on-site repulsion. The 
SDW fluctuations dominate at large U and near zero and 
half filling, where correlation effects between electrons 
are known to be important. The region in parameter 
space where triplet superconducting fluctuations domi- 
nate is separated from the region of SDW fluctuations by 
a metallic region. This constitutes a novelty, as previous 
attempts to draw a phase diagram for purely instan- 
taneous interactions predicted sharp transitions between 
the SDW (CDW) and superconducting regions. The re- 
tarded interaction mediated by the phonons is responsi- 
ble for this new phase. For a fixed U, we found that as 
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the phonon coupling is increased, superconducting fluc- 
tuations are most likely to occur near quarter filling, or 
away from half filling and zero filling where correlation 
effects are most important. A natural extension of our 
investigation would be to study other models of interact- 
ing electrons with nearest neighbor interaction, such as 
the extended Hubbard model, where, however, it is nec- 
essary to rely on numerical estimates for the Luttinger 
liquid parameters. 

We noted that by varying the electron filling factor, 
the Wentzel-Bardeen singularity can be approached for 
arbitrary values of the electron-phonon coupling. This 
singularity was previously understood to be the 

point where the coupling is so strong that the system ac- 
quires a negative compressibility. Nevertheless, for corre- 
lated electrons systems with on site repulsion such as the 
Hubbard model, the singularity is reached by decreas- 
ing the charge velocity associated with the particle-hole 
excitations. Near this point, SDW (CDW) correlations 
are totally suppressed, and the system is pushed through 
the metallic phase towards the superconducting phase. 
We finally note here that this behavior is reminiscent of 
a characteristic feature of high temperature supercon- 
ductivity materials: the increase of doping (and thus 
change of filling factor) changes the phase of the material 
from antiferromagnetic to metallic and finally to super- 
conducting (at low enough temperatures). Since there is 
some indication that the 2D layers of these new mate- 
rials can be described in terms of Luttinger liquids [Q , 
it is not unreasonable to conjecture that the mechanism 
discussed here could be of relevance for these materials. 
However, this is an open question which deserves further 
investigation. 

Next, we argued that our results can be applied to 
study the (in)stability of two electrons (or one electron- 
hole) modes in an electron wave guide, coupled with 
a short range density-density interaction. For specific 
Hamiltonians decribing the interaction between electrons 
within the chains, we can determine the threshold cou- 
pling where the chains become superconducting. Alter- 
natively, the transition to the superconducting regime 
could be observed by varying the electron density, expoit- 
ing once again the WB singularity. This could be probed 
by studying the periodicity of the persistent current of 
a two channel mesoscopic ring as a function of electron 
density. Such rings with few transverse channels can now 
be fabricated using GaAs/AlGaAs heterostructures with 
metallic gates, as in the remarkable experiment of Mailly 
et al. 1^^. The electron density can be varied my means 
of an electrostatic gate. The flux causes a twist in the 
boundary conditions, but docs not affect the bulk proper- 
ties, such as the correlation function exponents. In addi- 
tion, the mutual inductance of the two transverse modes 
can be neglected because its effects are proportional to 
the small ratio avp/ci^ with a ~ 137 the fine structure 
constant and c/ the velocity of light In the metal- 

lic regime, the power spectrum has a peak associated 
with the fiux quantum (f)Q — hc/e, whereas in the super- 



28 



conducting case, this peak should give place to a peak 
at the superconducting flux quantum 4'q/2. Finally, the 
study of the stability of two coupled transverse channels 
could in principle be generalized to an arbitrary number 
of channels, as it was recently shown that the correlation 
function exponents for the multicomponent Tomonaga- 
Luttinger model can be obtained formally [Q using con- 
formal invariance arguments. 

There are still challenging problems to be addressed. 
First, we have neglected throughout this paper the cou- 
pling between the fast oscillating part of the electron op- 
erator to the 2fci? phonon mode. This coupling is known 
to be at the origin of the Peierls instability, but is absent 
for the case of the coupled electron modes in the wave 
guide. For electron-phonon systems, it would be inter- 
esting to study whether the dramatic effects triggered by 
the Wentzel-Bardeen singularity can survive this addi- 
tional interaction. This question can be answered in the 
case where the electron bandwidth is comparable to the 
phonon cutoff, where the retardation effects associated 
with the 2kp phonon processes can be neglected, and this 
additional coupling amounts simply to a reduction of the 
on site repulsion parameter. The Wentzel-Bardeen sin- 
gularity survives the 2kp phonon processes in this case. 
In most physical systems, however, the phonon cutoff is 
smaller than the electron cutoff, and these retardation 
effects have to be taken into account . 

Second, the present treatment uses exact results for 
the isolated Hubbard model, and applies it to a system 
which is coupled to phonons. It would be desirable to 
have a self-consistent treatment, where the effects of the 
correlated electrons on the lattice are also described . 
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APPENDIX A: TIME EVOLUTION 



In order to calculate the commutator in Eq. (4.4), we 
need to find the explicit time dependence of the opera- 
tors Lpp, ipfj^ Hp and Hg.. The evolution in real time is 
determined by the following equation of motions {% = 1): 







(Ala) 






(Alb) 


dt^p 


= ^dl^p+gdld 
Kp 


(Ale) 


dtd 




(Aid) 




= Ceroid + gdlipp . 


(Ale) 
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Assuming for the above quantities a dependence in time 
and space of the form exjp(ikx — iujt) we obtain the dis- 
persion relations: 



~ vlk 



2 2 I 2 



(A2a) 
(A2b) 



for the charge and spin fields. The hybrid velocities v± 
are defined in Eq. (4.15). 

At t = 0, the fields can be expressed as a Fourier series: 



1/2 



np(x,0) = *5](^)'^V(4-a.,) 



1/2 



^^(-'O)-E(^) ^^^^^^^ 



e^""(6l. + 6-fc) 



2L 



1/2 



1/2 



(A3a) 
(A3b) 
(A3c) 
(A3d) 
(A3e) 
(A3f) 



where (bk) is an annihilation operator characteriz- 
ing a charge (spin) excitation at momentum fc, and Ck 
is a phonon annihilation operator. These expressions 
are implicitly regularized by the large momentum cut- 
off, exp{—e\k\L/2'!T), and furthermore the fc = mode is 
excluded in the sum. 



Eq. (A2a) indicates that the time dependence of the 
charge fields should be of the form: 



f 1 



\2L\k\ 



1/2 



Akx 



a,fi=± 



(A4a) 



2L 

a,f}=± 



(A4b) 



Similarly, the tim e dep endence of the spin fields is deter- 
mined from Eq. ( A2b ): 



1/2 



Akx 



(A5a) 
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E 

k 



\k[\ 
2L J 



1/2 



Akx 



D+e 



'ikUa- 1 



ikua- 1 



(A5b) 



The operators A±±^ B±±, C± and D± are determined 
by expressing the initial values ip,j{x,0), dt^v{x,0), 
d^(f^ {x, 0) , d f(p^{ x, 0), in terms of spatial derivatives via 
Eqs. (Ala)- (|A1e|), and similarly for IIj^. We then obtain: 



(3upKp 
1 — ^sgn 



(fc) 



^a-fc 1 + 



+iKpUpg{cC) 



Va 
-1/2 



sgn(fc) 



Cfe(-1 H 

Va 



/3c 



sgn(fc) 



(A6a) 



/3 



KpVaUp 



+.9(cC) 



-1/2 



-^«-^.2„2^ + 6>gn(fc)) 



Cfc -c K + - fl„)sgn(fc) 



+cU ( -c + — (Wp + - «-a)sgn(fc) 



Caik) = 
Da{k)=i 



1 - aK^sgnik) 



1 _ sgn(fc) 

2 " 2i^, 



1 + Q;ifo.sgn(A:) 



2 

sgn(/i;) 



2 2/C 



b-k 



(A6b) 

(A6c) 
(A6d) 



with a, P — ±. These operators satisfy the commutation 
relations: 



[Bai3{k),A^^{k')] = ~iSk,k'Sa,pS^^^ 

[Da{k),Cf3{k')] = -iSk,k'Sa.,l3 



(A7a) 
(A7b) 



These relations are used in turn to calculate the commu- 
tator s whi ch appear in the normal ordered products of 
Sec. IV A , with the substitution t ^ it for the thermo- 
dynamic quantities: 

[ep{x,T),^p{QM 



I f v^p-vt 



In 



1 — z{x — v+ir) 1 — z{x + v+ir) 
At: \ v'j^ -\- v'^ \_ \ — z{x — v+ir) \ — z{x + v+it) 

1 — z{x — V-ir) 1 — z{x + v^ir) 



9 9 

Up-vj 



In 



[0,(a;,r),(^,(O,O)] 



In 



1 — z{x — v^it) 1 — z{x + V-It) 

(A8a) 

1 — z{x — Ufjir) 1 — z{x + UcIt) 
1 — z{x — u^it) 1 — z{x -\- u^ir) 

(A8b) 
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Noting that [0^{x,t), (p^{0,0)] = [9^(0,0), ip^{-x, ~t)], 
we obtain the exphcit form for (;„. If in addition x ^ 0, 
we may put e = and use that z{x + viT) — 1 / z(x_-- vir) , 
which then leads to the second equahty in Eq.(4.5) after 
some simple manipulations. 



APPENDIX B: FOURIER SUMS 



In this Appendix we compute the Fourier sums which 
appear in Eq. (4.13) in the limit T 0. For that purpose 
it is convenient to rewrite the quantities to be summed 
over in terms of irreducible fractions expressed as func- 
tions of and k"^: 



2 9 

c — vt 



1 



1 



1 



Kpulk^Dpik) 



1^2 _|_ „2 ^2 / 

— — v'^ 

2 9 



(Bla) 



1 



1 



+ v^k'^ uj'^ + v^k'^ 



1 



;2 + v-^k'^ 



(Bib) 



where the velocities v± are defined in Eq. ( 4.15| ). 

With the above decompositions, the sums over Mat- 
subara frequencies in Eq. ( 4.13| ) yield, in the T — > 
limit: 



UpKp 



1 



ujsmiujT] 



-y 



1 ^ sin(a;T) 



2|/c|(u+ + W_) \ Vj^V- 

2\k\ 

C 



sgn(r) 



sgn(r)^e--l^^l 



V- \kT\ 



v'i ~ 



(B2a) 
(B2b) 

-v+\kT\ 

(B2c) 
(B2d) 



1 ^ cos(tJT) _ Kp f Up vl- 



Dp{k) 2\k\ Vw- vl~v 



9 9 

— V_ 



1 ^ cos(ijJt) K, 



E_Q-Ua\kT\ 



(B2e) 
(B2f) 



where sgn(r) = 1,0, —1, if r >,=,< 0. Finally, the sum- 
mations over momentum which occur in the calculation 
of the Green function exponent are: 



TT X ^ 1 

cos(fcx)e~^l'=^l 



Ine (B3a) 



ln|l - z(a; + ii;|r|)| (B3b) 
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TT s-^ sm{kx)e "l*^"^' 1 1 — z(x + zu|t|) 



— In ) -r^T , (B3c) 

2i l-z(x + iv\T\) ' ^ ^ 



k 2i l — z{x + iv\T\) 

where z{x) = exp(27r«a;/_L — e), and where again 
we have imphcitly used the large momentum cut-ofF, 
exp{—e\k\L/2n), and omitted the k = term. 
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FIG. 1. Green function exponent for the Hubbard model: 
a) at fixed filling factor, as a function of u = U/t for b/cup — 
(solid line), b/cup — 0.2 (dotted line), b/cUp = 0.4 (dashed 
line) and b/cup = 0.6 (dashed-dotted Une). b) for fixed U, as 
a function of filling factor and for the same values of b/cup as 
in Fig. la. 



FIG. 2. Phase diagram for U — O.Olt, for low filling factor 
n < 0.2, as a function of the phonon cou pling b/cUp, usin g 
the perturbative result of the inequalities (6.4a) and (6.4b). 
SDW, M, and SC refer to the spin density waves, metallic, 
and triplet superconducting regions. 



FIG. 3. Phase diagram of the Hubbard model coupled to 
phonons, as a function of filling factor and phonon coupling 
constant b/cUp-. a) U/t = 0.7. b) U/t = 0.3. c) U/t = 0.1 . 
We use the same labels as in Fig. 0. 
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FIG. 4. Phase diagram of the Hubbard model coupled to 
phonons at quarter filling (n = 1/2), as a function of U jt and 

6/cMp. 

FIG. 5. Plot of Up/Kp as a function of the filling factor n. 
Up is meeisured in units of vp{n — l)/2 = at. Prom top to 
bottom: U/t = 16,8,4,2. Note the abrupt change for small 
values of U as one approaches half filling. 

FIG. 6. Phase diagram for the extended Hubbard model 
in the limit U — > cxd, as a function of the nearest neighbor 
interaction V/2t and the phonon coupling constant h/cUp. 
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